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We introduce an interference measure which allows to quantify the amount of interference present 
in any physical process that maps an initial density matrix to a final density matrix. In particular, 
the interference measure enables one to monitor the amount of interference generated in each step 
of a quantum algorithm. We show that a Hadamard gate acting on a single qubit is a basic 
building block for interference generation and realizes one bit of interference, an "i-bit". We use 
the interference measure to quantify interference for various examples, including Grover's search 
algorithm and Shor's factorization algorithm. We distinguish between "potentially available" and 
"actually used" interference, and show that for both algorithms the potentially available interference 
is exponentially large. However, the amount of interference actually used in Grover's algorithm is 
only about 3 i-bits and asymptotically independent of the number of qubits, while Shor's algorithm 
indeed uses an exponential amount of interference. 

PACS numbers; 03.67.-a, 03.67.Lx 

I. INTRODUCTION 

Entanglement and interference are believed to be key ingredients that distinguish quantum from classical information 
processing jj. Indeed it has been shown that large amounts of entanglement must necessarily be generated in a 
quantum algorithm that offers an exponential speed-up over classical computation Q. Tremendous effort has been 
spent to develop methods to detect and quantify entanglement in a given quantum state (see 0, Q for recent reviews) . 
On the other hand very little has been done to quantify "interference" . It seems therefore well worthwhile to analyze 
the role of interference in quantum algorithms in more detail. 

While entanglement plays an important role in many quantum information tasks, like for example quantum tele- 
portation , quantum key distribution y , or super-dense coding J3, large amounts of entanglement may not be 
the only requirement to get a speed up with a quantum algorithm j^, l3l ■ Classical analogues of entanglement exist 
for example in the context of the propagation of classical phase space density through Liouville's equation [lol |. 
Even without talking about dynamics, a formal analogy of pure state quantum entanglement can be easily defined in 
the classical domain, if we replace wave functions in the definition of quantum entanglement by probability distribu- 
tions. Arbitrary large amounts of "classical entanglement" of many-particle probability distributions (corresponding 
to classical correlations) can thus be created by a classical stochastic computer. However, only linear combinations of 
probability distributions with real positive coefficients are allowed, and this makes it impossible to efficiently realize 
Fourier-transforms with high contrast, a decisive ingredient in many quantum algorithms (see e.g. |lllll2l[T3j ). In gen- 
eral it seems that interference between many computational paths plays an important role in quantum computation. 
It is unclear, however, how much interference, if any, is needed for a given quantum speed-up. 

As a first step in the direction of answering this question, we introduce here a general interference measure that 
applies to any physical process which maps an initial density matrix to a final density matrix. It can be used in 
particular for any quantum algorithm, measurement processes, or classical communication. We use the interference 
measure to quantify interference in various steps of the well-known quantum algorithms of Shor and Grover [l^ 
as well as other physical examples, including photons passing through a beam-splitter or through a Mach-Zehnder 
interferometer, decoherence through bit-flip and phase errors, and quantum teleportation. 

II. THE ESSENCE OF INTERFERENCE 

The familiar example of a double slit experiment, where waves with different wave- vectors are superimposed at the 
slits, propagate and interfere to generate an interference pattern, i.e. a certain probability distribution on the screen 
in a position basis, contains the basics of any interference experiment: 

1. Interference needs coherence. An interference measure has to distinguish between coherent and incoherent 
propagation. This means first of all that interference is a property of a propagator of states, not of a state itself, 
in contrast to entanglement. 

2. Coherent propagation alone does not amount to interference, however. For example a quantum gate which 
implements just the unitary identity operation is completely coherent (a pure state remains a pure state). 
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whereas there is no interference at all: the probability for each final basis state depends only on the initial 
amplitude of the same basis state. Similarly, there is no interference for a quantum gate that just permutes 
in-coming amplitudes. Evidently, at least two states have to be coherently superposed. Generally speaking, the 
amount of interference should be linked to how many different initial state amplitudes contribute to each final 
state amplitude, and to what extent. An interference measure should be maximized if each basis state as input 
state produces an "equipartitioncd" output state, i.e. a state with the same absolute probability amplitude for 
each basis state. 

3. Interference is basis dependent. Indeed, no interference pattern at all, but a probability distribution just 
reflecting the initial amplitudes would be observed in the double slit experiment described above, if the output 
were observed in the momentum basis. Also any unitary evolution is given by a unitary matrix U defined in 
a specific basis (for example the computational basis for a quantum algorithm). This matrix can always be 
diagonalized, and in the new basis the propagator then just corresponds to a multiplication with a phase factor 
for each input basis state, which does clearly not amount to interference. 



III. THE INTERFERENCE MEASURE 



A. General formulation 



We start by defining a measure of coherence for a general propagator P of a density matrix, p' = Pp. where P is a 
super-operator specified in the computational basis, such that, written in components, 

Pij =^Pi]MPkl ■ (1) 

k,l 

Further below we will adapt this formalism to the "operator sum formalism" more familiar in the context of quantum 
computation [l5j |. 

Consider a black box which maps an initial density matrix to a final probability distribution, as is normally the 
case for a quantum computer if the read-out process is considered part of the algorithm. As long as one regards a 
unique input state (pure or mixed) , one cannot distinguish between coherent and incoherent propagation. Indeed, the 
black box might forget about any initial phase information, keep just the initial probabilities and then use a stochastic 
matrix which maps them to the probability distribution of the desired final state. The final probabilities are then 
invariant under arbitrary phase rotations of the initial amplitudes. To distinguish between coherent and incoherent 
propagation, an interference measure must therefore quantify the dependence of the final probabilities on the initial 
phases — a strategy which is in fact often employed experimentally to show coherence. 

Let us therefore start with a pure initial state, p = \il'){ip\, with \ip) — The amplitudes Oj with phases 

(fij, ttj = |aj |e"^J , lead to final probabilities 

P'.^lpul - (2) 

The dependence of each individual probability on the initial phases is thus given by 

1^ = i^|afcan(P.Mfce'(^'~'^'')-c.c.) . (3) 

We define the real phase sensitivity matrix S with matrix elements Sn = dp'^/dtpi, and the positive semi-definite 
matrix S'S'"^, with SS'^ = iff dp^/dipi = for all «, ^ = 1, . . . , TV. A coherence measure C(P) can be obtained by 
taking the trace of SS'^ and averaging it over all initial phases, 

C(^) = / d^^...d'fNtr{SS'^) (4) 

r27r i_ i_ T„ ,.2ir 



(27r)^ 

.»,fc,;,m L Jo K^^) Jo l^^j 



= 2EE|P,,,ifcafeaip. (5) 

i,k ly^k 
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This measure still depends on the amplitudes of the initial state. However, coherence should be measured for an input 
state with amplitudes on all basis states. We therefore chose a "democratic" input state with |ai| = 1/\/N for all 
i = 1, . . . , iV. We multiply with an additional prefactor N"^ /2 and define the coherence measure 



X(P) = ^|P, 



ii,kl I 



■i.k.l 



i,k 



kk\ 



(6) 



The quantity T{P) obviously has the desired property to be zero in the case of a classical stochastic propagation, 
that is for Pij^u = MikSijSki, where Mik is a stochastic matrix which propagates classical probabilities from initial 
values pk to final values p'^, p'^ = J2k ^ikPk, and Sik denotes the Kronecker-delta. From the definition it is clear that 
X(P) is non-negative. In the case where all eigenvalues of P are smaller than one, as for example for dissipative 
quantum maps an upper bound is given by I(P) < N'^, as can be seen by using the inequality UPV"!! ^ IIV'II with 
V'fez = Sk,ko^i,io for all 1 < fcoj^o ^ N. This bound is probably not optimal, and, as will be seen, can be improved in 
the case of unitary propagation. 

The propagation of mixed states is in quantum information theory generally formulated within the operator sum 
formalism ^ set of operators {£';} acts on p according to p' = EipEj — Pp. The i?/'s are known as Kraus 
operators |T9| and obey 'J2k -^l-^k = 1 for trace-preserving operations. The connection to the propagator P is given 
by Pij^krn ~ '^i{Ei)ik{Ei)jjm ^nd we can reformulate the interference measure 10 as 



I jirn 



I ik\ 



(7) 



We will now show that I{P) is in fact a good measure of interference, as it also measures the amount of equipartition 
in the case of unitary propagation. 



B. Unitary propagation 

Coherence is, by definition, perfect, if all pure incoming states remain pure during propagation, i.e. also the final 
state can be written as p' = |$)($|, with a state |<i>) related to \ip) by a unitary transformation, 1$) = U\tp). In this 
case P, which we shall denote by P{U), has matrix elements Pu.ki — UikU*i. Profiting further from the fact that 
X^ffe 1=1 \UiiUik\^ = TV due to the unitarity of U , we find in this case 

Apm = ^^-Y.\u^k\'^ ■ (8) 

In the unitary case the coherence measure (jHJ has the property that < I{P{U)) < N — \. The right hand side 
of this inequality is easily verified using the Cauchy-Schwarz inequality applied to all vectors . . . , jf/iArp) 

(i = 1, . . . , iV), and (l/ViV, . . • , l/-\//V), and is in fact the optimal bound, as it is reached for \Uik\ — 1/VN for all 
i, k. The left hand side follows from the positivity of SS'^ , or exphcitly from ^ ■ \Uik\^ < J^i k Wik]"^ ~ J2k ^ — 
where we have used that < \Uik\'^ < 1 due to J^i \ Uik\'^ = 1. 

The term is nothing but the inverse participation ratio (IPR) of a column k of the unitary matrix 

summed over all columns. The IPR is a well-known measure of "equipartition" of a wave-function, used extensively 
in solid state physics in order to measure localization jl7i |. A column of U with an amplitude on a single basis state, 
i.e. Uik = 5ii^, for column k and some index ik gives |CAfc|^ = 1. If all columns have an entry on just a single basis 
state we get therefore ^^k Wik\'^ — ^ ■ 

Thus, if Vik = Sip.(k), where Pi{k) is a permutation, we have I{P{U)) = 0. This reflects just the fact that this 
kind of coherence is useless — the final probability distribution does not depend at all on the initial phases and for 
all possible input states, the same output could be obtained with a propagation of probabilities only. 

On the other hand, perfectly "equipartitioned" output states for each computational basis state used as input, 
\Uik\ — I/'n/ZV for all z,fc, give ^j^^, iC/ifcj^ = 1. Therefore, I measures for unitary propagation also the amount of 
equipartition, where X{P{U)) varies between for a mere permutation of all input states and A^ — 1 if all output states 
corresponding to computational basis input states are perfectly equipartitioned. We therefore define T{P) in eq. 0) 
as the measure of interference in the propagator P. Note that this definition is very general, and might be applied to 
any physical situation where (possibly mixed) input states are propagated to (possibly mixed) output states. 
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C. Properties of the interference measure 

1. Invariances 

As eq.® contains a double and triple sum over all computational basis states, it is obvious that I is invariant 
under a permutation of the computational basis states. This goes hand in hand with the observation that a pure 
permutation of the computational basis states does not generate any interference. It is also obvious that I is invariant 
under arbitrary phase changes of any matrix element Pij^ki- For unitary propagation, the inverse propagation 
leads to the same amount of interference, l{P{W j) = X{P[U)). 

2. Equipartition in a Sub-Space 

Consider a unitary matrix with \Uik\ = I/a/M for 1 < i, fc < M for some integer M < N, and \Uik\ ~ 6ik for i > M 
or k > M . Straightforward evaluation of T{P{U)) from eq.ljHJ gives I{P{U)) = il/ — 1. Thus, I{P{U)) increases 
indeed linearly with the number of coherently superposed states. 

3. Adding auxiliary qubits 

Another situation, encountered often in quantum computing, is the addition of auxiliary qubits, or in general a 
Hilbert space of dimension M connected by a tensor product to the original Hilbert space. As long as one acts on the 
auxiliary qubits only with the identity operation, their only effect is to increase the number of coherently superposed 
states by a factor M. One should therefore expect the amount of interference of C/ = U (S) 1m to be larger by a factor 
M compared to X{P{U)). This is indeed the case. To see this, define the matrix elements of JJ as Unmki = Unk^mi = 
UnkSmi- Then we have from eq.©, I{P{U)) = NM - En,m,k,i \Unk\^Smi = M{N - J2n,k Pnkl^) = MI{P{U)). The 
corresponding calculation can be done for the general non-unitary propagation of a density matrix p — p"^ ® , such 
that p' = Pp= (Pp^) «) p^, and the result is the same: J(P) = MI(P). 

4- The "i-bit" and the Hadamard gate 

The linear dependence of the interference measure on the number of coherently superposed states makes it possible 
to define a unit of interference. We can define the number nj of "interference bits" (or "i-bits" for short) that measures 
the (logarithmic) amount of interference in a propagator P as 

n, = log2(X(P) + l). (9) 

As an immediate consequence we find that a Hadamard gate, Hij — for i,j £ 0, 1, generates one bit of 

interference as T{P{H)) = 1, and this is the maximum possible amount of interference for the unitary propagation 
of a single qubit. Moreover, it is easy to show that for a tensor product of n Hadamard gates acting in parallel on n 
qubits, also called the Walsh-Hadamard transform, W — H®H®...®H, one obtains n i-bits, X{P{W)) = 2" — 1, 
as the amplitude of each matrix element of the tensor product has an absolute value 1/2"/^. One may thus consider 
the Hadamard gate as a "currency" of interference in a quantum algorithm, much as a singlet measures the amount 
of entanglement in a bipartite quantum state. Note, however, that not each Hadamard gate in a quantum algorithm 
adds an i-bit to the total amount of interference. If p Hadamard gates act on p different qubits out of a total number 
of n qubits, the interference is I = 2" — . 

D. Potentially available interference versus actually used interference 

Another reason why a Hadamard gate does not necessarily add an i-bit of interference to a quantum algorithm lies 
in the fact that the amount of interference added in a given step of the algorithm depends on the transformation built 
in all previous steps. This is easily seen from the example of two Hadamard gates in series, acting on a single qubit, 
= I2, such that T{P{H^)) = 0. We will call "accumulated interference" the total interference T{Pm ■ ■ ■ Pi) of a 
sequence of transformations P,„ • • . -Pi , which is in general very different from the sum of interferences of each step 
Pi. In principle one can calculate the accumulated interference for an arbitrary part of any quantum algorithm, but 
most interesting is in general the accumulated interference of the entire algorithm. 
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It turns out that many quantum algorithms generate a lot of interference right at the beginning, including Shor's 
and Grover's algorithms (see below), as they start out with the Walsh-Hadamard transform on many qubits. In the 
end the desired information (e.g. the period of a function or the label of a searched state) has to be extracted, which 
is done by a reduction of the accumulated interference, as the probability distribution gets concentrated on only very 
few computational states. The Fourier transform, which taken by itself would add even more interference, performs 
this task in Shor's algorithm. 

However, not all quantum algorithms actually use all the interference generated. A tremendous "waste" of interfer- 
ence can arise when only a single column of the unitary matrix corresponding to a single initial state (e.g. |0 . . . 0)) 
is used, whereas it is completely irrelevant what happens in the other columns. One might be tempted to calculate 
the interference just corresponding to that initial state. Formally this can be done by including a projection onto 
the state |0. . .0) in the algorithm, but it is easily shown that then the accumulated interference of this projection 
combined with whatever unitary operation follows vanishes. All coherence is destroyed in the sense that the final 
probability distribution is independent of the single remaining initial phase. An alternative approach is to calculate 
the accumulated interference for the remaining algorithm U after the initial Walsh-Hadamard transform, i.e. once 
amplitudes from all computational states are available in the superposition. This amounts to changing the initial 
state. Then all columns of U do contribute, and as a consequence X{P{U)) is a more realistic measure of the interfer- 
ence actually used. The interference measures corresponding to U and U therefore give complementary information 
(total "potentially available" interference in the algorithm and "actually used" interference, respectively) and will be 
calculated below for Grover's and Shor's algorithms. 



IV. APPLICATIONS 
A. Beam splitter 



A beam splitter transforms two photon modes with annihilation (creation) operators a,b (a^,6^), respectively, 
according to the unitary transformation Ubs = exp(0(a^6 — ab^)) |l5j . The action of Ubs on a state \nm) with 
m photons in mode a and n photons in mode b can be easily derived with the help of the relations t/BS^C^Bs — 
a cos 9 + bsinO and C/rs^C^bs — bcosd — asinO. We exploit the fact that the total number N of photons in the two 
modes is conserved to write the matrix elements of Ubs as 



{Uif)„n = {i N ~ i\UBs\m N - m) 

j ., . ~! Min{N—i,m) / \ / ^.j \ 

l^m^ ^ (m\( N-m ^(_iy(,o,^)™+^-.-2.(^i^^).-™+2Z, (10) 



V ^ ' i=Max(m-i,0) ^ ^ ^ 

For example, the dual-rail representation jl^] with logical states |01) and |10) leads to the 2x2 rotation matrix 

UW=(^0S9 -sine\ 

ysm9 cosd J ' ^ ' 

with the amount of interference J{P{UBg)) = 2(1 — cos'' 6* — sin'' 6*). Thus a maximum interference of 1 i-bit can 
be achieved for 6 = 7r/4, in which case the beam splitter realizes indeed just a Hadamard gate, up to phase shifts. 

The function 2{P{uj^\6))) appears to be periodic with period 7r/2 for all N. The maximum amount of interference 
increases approximatively linearly with N and is reached for sufficiently large iV at = 7r/4. With increasing TV, the 
maximum is reached more and more rapidly as function of 0, and X remains almost flat in a broader and broader 
interval around the maximum (see figure 



B. Mach-Zehnder interferometer 



A Mach-Zehnder interferometer consists of two beam splitters in series, with an additional phase shifter in one of the 
two arms [T5l |. The second beam splitter is inverted relative to the first one, such that the total unitary propagation 
for N photons is given by 



(12) 
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FIG. 1; (Color online) Amount of interference 1 generated by a beam splitter as function of the angle Q and A'' = 1, ... 20 
photons (top). The amount of entanglement at Q = 7r/4 increases approximately linearly with A*" (bottom). The full line is 
T = N. 

The phase shifter acts on a m photon state in mode a as \J^\nm) — e"^™|nm). The interference generated by the Mach- 

Zehnder interferometer is periodic in Q and with periods 7r/2 and 27r, respectively (Figure [21). Again, l{P(ul^^)) 
becomes more and more flat around the maximum with increasing TV. The settings 6 = tt/2, or 9 = 7r/4 or 37r/4 and 
4> = IT, lead to a minimum of T in the range of N investigated (N — 1, . . . , 20). Finally, = never leads to any 
interference, as in this case the interferometer just performs the identity operation. 

C. Decoherence: bit flip errors and phase errors 



The error operators for a bit flip error in a single qubit are given by Eq — ^/pls, Ei = ^1 — pX, where X is the 
Pauli ax matrix and p the probability for no error to occur. Similarly, the error operators for a phase error in a 
single qubit read Fq = ^yp'i-2, Fi = V^~-pZ, where Z is the Pauli az matrix. One easily shows using eq.|(Zj) that the 
amount of interference vanishes for both types of errors, as it should of course be, due to the fact that these are purely 
decohering processes. More interesting is the situation where we first apply a Hadamard gate and then the errors. 
That is, the Kraus operators are now EqH,EiH for a bit flip error after the Hadamard gate, and FqH,FiH for a 
phase error after the Hadamard gate. In the first case, we obtain T — {1 — 2p)^. As it should be, the interference is 
completely conserved, X = 1, for p = or p — 1, i.e. when bit flips either never occur or do occur with certainty. Even 
in the latter case the interference is perfect, as the bit flip just corresponds to a unitary permutation. On the other 
hand, in the case of a phase flip after a Hadamard gate, interference is always conserved, i.e. 2=1 independent of 
p. While this might be surprising at first, it is in fact a remarkable property of I to detect that the phase error does 
not modify at all an interference pattern obtained by coherently superposing |0) and |1). To see this, we start with 
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FIG. 2: (Color online) Accumulated interference generated in a Mach-Zehnder interferometer as function of the angles 9 and 
cj> for 1 photon (top) and 20 photons (bottom). Note the different scale. 

the initial state (e''^|0) + |l))/\/2. After the Hadamard gate this state becomes (e''^(|0) + |1)) + (|0) - |l))/2. With 
probabiUty p the state remains the same, while with probability 1 — p a, phase error flips the sign in front of the |1) 
states. One obtains thus the final density matrix p/ = cos2((y5/2)|0)(0|+sin^(</j/2)|l)(l| + (l/2-p)isin¥3(|l)(0|-|0)(l|). 
Clearly, the probabilities to find or 1 in the final state are unaffected by the phase error, and one gets a perfect 
interference picture for any value of p, as the phase errors from the bra and the ket cancel. Note that by just looking at 
the off-diagonal matrix elements of p/ (as is often done to estimate the amount of coherence in a state) one concludes 
that complete decoherence and the reduction to a classical mixture occurs at p = 1/2, whereas for p = Oorp=la 
completely coherent final state is retained. This is, however, irrelevant for the success of the interference experiment, 
and correctly detected by our interference measure. 

More generally, one might consider the action of the phase error as a measurement of the interference pattern. In- 
deed, Z is diagonal in the computational basis, which thus constitutes the "pointer basis" of Z 18] . The measurement 
reveals the interference pattern, but does not destroy it. It is therefore reasonable to demand that an interference mea- 
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FIG. 3: (Color online) Accumulated interference during the quantum teleportation protocol |3- On the t-axis is the step 
number (see text). The interference is plotted for the entire propagator up to and including step number t. Interference only 
arises during the application of the two Hadamard gates and is conserved during all other steps, including the measurement 
(step 5). A maximum amount of interference of about 2.58 i-bits (X = 6) is generated. 



sure be conserved during the measurement process. The fact that interference can be conserved during measurements 
is crucial also for quantum teleportation. 

If a phase error occurs after the phase shifter in the Mach-Zehnder interferometer, the interference is reduced. For 

the example of one photon the Kraus operators are given by Go — VPf^BS ^p^bs'^ ^^'^ ~ Vl ~ P^^bs'^^p^bs ^- 
The general expression for the interference is somewhat cumbersome, but simphfies for beam sphtters with 9 — n/A 
to T = (sin(y9(l — 2p))^. So the interference in the Mach-Zehnder interferometer is clearly completely destroyed for 
p = 1/2, as is to be expected. 



D. Quantum Teleportation 

Alice can tclcport an unknown quantum state of a qubit 1 to Bob, by doing the following 5j: She first prepares 
two auxihary qubits (2,3) in the Bell state (|00) + |11))/V2, and sends qubit 3 to Bob. We will consider the Bell state 
preparation as part of the protocol, as it will become clear that already here a certain amount of interference is used 
|23l |. Starting from state |00) of qubits (2,3), Alice can prepare the Bell state by applying H to qubit 3 {H3, step 1) 
and then a CNOT with 3 as control and 2 as target (C32) to qubits (2,3) (step 2). Next Alice applies a CNOT with 1 
as control and 2 as target (C12) to qubits (1,2) (step 3) and a Hadamard to qubit 1 {Hi, step 4). She then measures 
qubits (1,2) in the computational basis (step 5), and sends the result mi, 7712 (where to^ = or 1) to Bob through a 
classical channel. Bob performs the unitary operation j^ni2 (gtep 6). The total propagation of the density matrix 
thus reads 

p'^ J2 Z^X'E,jHiC\2C32H3p{Z^X'E,jHiCi2C32H3y , (13) 

where the Eij are projectors, Eij = 1 (g) \ij){ij\, and the qubits are numbered from right to left. Figure |3| shows how 
the amount of interference used in this protocol evolves step by step. Obviously already the creation of the initial 
entangled state uses interference, through the application of the Hadamard gate in qubit 3 (step 1). Qubits 1 and 2 
are left alone, therefore the total amount of interference is I{P{H (g) I4)) = 4. In fact, interference is only generated 
by Hadamard gates in this algorithm, and the largest amount of interference reached is X = 6, or almost three i-bits. 
Remarkably, the interference is not destroyed in the measurement process. 



E. Shor's algorithm 

Shor's algorithm (llj enables one to factor a large integer number R into primes using a polynomial number of 
operations. The initial state is prepared in two registers of size 2L and L, where L = [log2 R] -\- 1 ([.] denotes the 
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FIG. 4: (Color online) Accumulated interference generated during Shor's algorithm for factorization of i? = 15 (L = [logj R] + 
1—4 and in total n = 3L = 12 qubits). On the t-axis is the step number (see text). The interference is plotted for the entire 
propagator up to and including step number t. Values of a are a = 13 (circles, full red line), a — 7 (squares, dashed green line), 
a = 11 (diamonds, full blue line), a = 8 (triangles up, full black line), a = 14 (triangles down, full cyan line), a = 4 (crosses, full 
orange line), a = 2 (stars, dashed violet line). Data for different values of a differ only for the last two time steps. Data from 
a = 13 and a — 7 are the same, as are data from a = 8 and a — 2. The horizontal dashed red line is the maximum possible 
value of interference for an untouched second register, X = 2" — 2^ = 4080. The inset shows the same curves for the last steps 
on a different scale. Lines are there to guide the eye only. Massive interference X = 2" — 2^ (or almost n i-bits) is generated 
during the application of the Walsh-Hadamard gate (part 1, first eight points). Interference is unchanged in part 2 (next eight 
points), and decreases during the last part (QFT) (last eight points). 

integer value), with a total number of qubits n = 3L. The algorithm can be decomposed into three parts. First, the 
state |0 . . . 0)|0 . . . 01) is transformed into N^^^"^ J^t^o^ \^)\^) (where N — 2^^ equals the dimension of the Hilbert 
space of 2L qubits) by the use of 2L Hadamard gates applied to each qubit separately. As mentioned in section 
nil C 41 this part generates the amount of interference, T = 2" — 2^. This value of X corresponds to the maximum 
value of I for transformations of the first register alone. On the other hand, no entanglement is created since the 
state remains factorizable. In a second part, the state is transformed in 0{n^) operations into N~^^^ X^tlo^ 
where f{x) = (mod R) . This part can be viewed as a permutation matrix: each line and each column has only 
one nonzero entry. Therefore the interference measure for this transformation alone gives zero. On the other hand, it 
creates entanglement. The third part consists of a Quantum Fourier Transform (QFT) on the first register only which 
allows one to find the period of the function /. The corresponding operator on the first register can be written as a 
matrix whose entries are all of absolute value 1/\/N. As such, the QFT alone maximizes the interference measure 
© on the first register, corresponding to X = 2" — 2^ , in the same way as the Walsh-Hadamard transform. In 
contrast with the latter, the QFT contains two-qubit gates and creates some entanglement. We note that nevertheless 
numerical simulations for small number of qubits seem to indicate that most of the entanglement is created during 
the second phase where no interference is built up or used .,19,] . 

It is interesting to note that various versions of the factorization algorithm have been proposed which aim at 
minimizing the number of gates and/or qubits poll2l|. They often use quantum Fourier transforms to gain efficiency 

in the process \x)\0 . . — > N~^/'^Y^^=I'o \^)\f{^)) ■• foi" example using instances of the Schonhage- 

Strassen algorithm which is more efhcient than usual multiplication for very large numbers. Thus these methods 
generate some interference in the second part of the algorithm, where usually only entanglement is produced. In some 
sense, this "acceleration" of the factorization algorithm is made by trading interference and entanglement. 

Figure 21 shows the accumulated interference for the factorization of i? = 15. Shor's algorithm requires a to be 
chosen coprime to R. We therefore show data for the seven possible values a = 2, 4, 7, 8, 11, 13, 14. The initial Walsh- 
Hadamard transform corresponds to the first eight value of time (the interference measure is calculated after each 
Hadamard gate). As explained in IIII C 4l after k Hadamard gates the value of X is X = 2" — 2^~^ . The whole 
transform generates an exponential amount of interference X = 2" — 2^ = 4080, which corresponds to the maximum 
value possible for a transformation of the first register alone. 
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FIG. 5: Accumulated interference up to and including step number t generated during Shor's algorithm, excluding the initial 
Walsh-Hadamard transform for factorization of i? = 15 (n = 12 qubits). The step number on the t-axis is the same as in figure 
l^but shifted by 8. Horizontal dashed red line is maximum possible value of interference for transformations of the first register 
alone I = 4080. The accumulated interference is changed only during the final quantum Fourier transform, and does not 
depend on the value of a. The actually used interference shown here is comparable with the potentially available interference 
shown in figure |31 

The next eight time values correspond to the construction of f{x) = (mod R) in the second register. Each time 
value corresponds to multiplication modulo R of the second register by , controlled by the value of the qubit of 
the first register. Each operation corresponds to a permutation in the computational basis, and does not affect the 
accumulated interference. 

The last part (QFT) corresponds to the last eight time values in figureSJ Each time step corresponds to a Hadamard 
gate followed by controlled phase transformations, which together build the QFT. Although the QFT taken alone 
generates interference, its net effect is to decrease the total accumulated interference by a small amount which depends 
on a. This is actually reasonable, since this last part concentrates probability on a certain number of states which 
depend on the value of a and R. It is interesting to note that the first six steps give the same values of I for all 
values of a. Only in the last two steps (corresponding to the least significant bits in the first register) does one see a 
branching. The first one distinguishes between the two values of the period of f{x) = (mod R) (the period is 2 
for a = 4, 11, 14 and 4 for a = 2,7, 8, 13), and the last operation distinguishes further between values of a. Data from 
a = 13 and a = 7 are the same, as are data from a = 8 and a = 2, which is consistent with the fact that for these 
values of a the set of values of (mod R) is the same. Other values of a correspond to different interference even 
though the period of f{x) may be the same since the total interference depends not only on the period of f{x) but 
also on the period of yf{x). The total interference which was exponentially large after the first two parts presumably 
remains so after the QFT, since for i? = 15 it goes down only by up to 7%. Note that it can go down and up in the 
last operations of the QFT. 

A more realistic picture of the accumulated interference "actually used" is obtained by calculating the interfer- 
ence measure for the operator obtained by omitting the initial Walsh-Hadamard transform (see section III D). This 
corresponds to the same algorithm but with initial state -^""'^^^ X]t=o^ This is shown in figure |31 where the 

time steps thus correspond to the last sixteen time steps in figure^ The interference measure is zero after the first 
part (construction of f{x) = (mod R) in the second register), which is in accordance with the fact that all eight 
steps corresponding to this part can be understood as permutations. On the contrary, the interference measure grows 
exponentially in the last (QFT) part, and reaches the maximum value 2" — 2^ (for transformations of the first register 
only). Thus the actually used interference is clearly exponential for the factorization algorithm. 
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F. Grover's algorithm 

Grover's algorithm Ug |0| can find a single marked item in an unstructured database of N items in 0{\/N) 
quantum operations, to be compared with 0{N) operations for the classical algorithm. The algorithm starts on a 
system of n qubits (Hilbert space of dimension N = 2"-) with the Walsh-Hadamard transform W, thus by building 
an uniform superposition of the basis states N^^/"^ Yl^l^o 1^)- Then the algorithm iterates k times the same operator 
U = WR2 WR1 w ith an optimal value k = [7r/(46l)] and sin^ 9 = l/N i.e. Ug = {WR2WRifW. As mentioned 
in section IlII C 41 W generates, taken for itself, the maximum possible amount of interference {n i-bits). The oracle 
Ri multiplies the amplitude of the marked item (the one searched by the algorithm) with a factor (—1), and keeps 
the other amplitudes unchanged. The corresponding N x N matrix is diagonal, making evident that its interference 
measure is zero. The operator R2 multiplies the amplitude of the state |0 . . . 0) with a factor (—1), keeping the others 
unchanged. By the same argument as for it generates no interference. Therefore all interference in the algorithm 
is generated by the Walsh-Hadamard transforms W. It is interesting to note that both operators i?i and R2 can 
be understood as multicontroUed gates which create entanglement, whereas W cannot create entanglement as it is 
composed of one-qubit operations. So Grover's algorithm alternates entanglement creation and interference creation 
during its evolution. We note that the interference is independent of the label of the searched item, as the interference 
measure X is invariant under permutation of the computational basis states. 

FigureElshows the evolution of the accumulated interference during Grover's algorithm in n qubits. In this figure, to 
avoid repetition, we took the Walsh-Hadamard transform as one single time step, contrary to the preceding section. 
Thus step 1 corresponds to the first application of the Walsh-Hadamard transform, which generates the massive 
interference I — 2" — 1; step 2 is the first application of the oracle Ri, step 3 the first application of the diffusion 
matrix D = WR2W, step 4 the second instance of the oracle R2, and so on up to a total number of 2fc -I- 1 steps. As 
expected, the oracle does not change the amount of interference, but the diffusion matrix reduces the accumulated 
interference in each step. This is crucial for the functioning of the algorithm: the probability flow is engineered in 
such a way that all probability gets concentrated on the computational basis state corresponding to the searched item. 
Therefore, the accumulated interference has to decrease, as the equipartition is decreased. 

As mentioned before, most of the interference built up during the application of the W gate is "wasted" , as U only 
acts on the initial state |0 . . . 0). Indeed, no particular state is selected by any other column of U. The accumulated 
interference therefore reduces only very slightly, from X = 2" — 1 down to X ~ 2" — 2. The situation is somewhat 
different for Shor's algorithm where the concentration of probability concerns many columns. 

As in the case of Shor's algorithm, it is also interesting to compute the actually used interference, i.e. the interference 
measure for Ug obtained by omitting the initial W transform, Ug = {WR2WR1)'' (see section HI D). Figure [7| 
shows that this number is much smaller than the potentially available interference. Moreover, it is asymptotically 
independent of the number of qubits. Indeed the interference for the first application of D is easily calculated 
analytically using the matrix elements = 2/N — Sij in the computational basis, and gives T{P{D)) = 8 — 
2A/N + 0{1/N'^). Thus, Grover's algorithm actually uses only about three i-bits, whatever the number of qubits on 
which it runs! After the first application of D, the interference shows damped oscillations with each application of 
the diffusion matrix, whereas the oracle Ri leaves the interference unchanged. Note that the oscillations would be 
undamped without the intermediate Ri steps, as {WR2W)^ = 1. 

Thus for the Grover algorithm the maximum actually used interference after the creation of the initial equiparti- 
tioned state is small and basically independent of the number of qubits. This is in sharp contrast with the results for 
the Shor algorithm (preceding section) where the actually used interference grows exponentially with the number of 
qubits. One may speculate if this is the decisive difference that leads to exponential acceleration of Shor's algorithm 
versus the acceleration for Grover's algorithm compared to the corresponding best classical algorithms known. 

V. CONCLUSIONS 

We have introduced a general measure of interference, which allows to quantify interference in any physical situation 
which involves the propagation of a density matrix. We have defined a logarithmic unit of interference, the "i-bit" 
and have quantified how much interference arises in each step of the two best known quantum algorithms as well 
as in other physical examples, including beam splitter, Mach-Zehnder interferometer, and quantum teleportation. A 
beam splitter and a Mach-Zehnder interferometer generate an amount of interference proportional to the number of 
photons, and quantum teleportation of one qubit needs about 2.58 i-bits. Both Shor's and Grover's algorithm build 
up an exponential amount of "potentially available" interference. However, Grover's algorithm actually uses only 
about three i-bits, asymptotically independent of the number of qubits on which it runs, whereas Shor's algorithm 
indeed uses an exponential amount of interference. It is therefore tempting to attribute the respective exponential 
versus VTV acceleration of these two algorithms to the amount of interference actually used, but more work will be 
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FIG. 6: (Color online) Accumulated interference in Grover's algorithm Ua 13 for n — 8 qubits. On the t-axis is the step 
number (see text). The interference is plotted for the entire propagator up to and including step number t. Massive maximum 
interference Tmax = 2" — 1 = 255 (or n = 8 i-bits) is generated during the application of the Walsh-Hadamard gate (step 1) 
and the accumulated interference decreases subsequently during the iteration of oracle and diffusion to the value X ~ 2" — 2, 
at which point virtually all probability corresponding to the initial state |0 . . . 0) is concentrated on the searched state. 
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FIG. 7; (Color online) Accumulated interference up to and including step number t generated during Grover's algorithm 
Og, excluding the initial Walsh-Hadamard transform [l^ for n — 7 qubits. The step number on the t-axis is shifted by 
one compared to figure |S] i.e. step 1 is now the first application of the oracle Ri, which does not lead to interference. The 
accumulated interference is changed only during the diffusion steps D — WR2W. The "actually used" interference shown here 
is much smaller than the "potentially available" interference shown in figure [5] and the maximum value X ~ 8 after the first 
application of D is asymptotically independent of the number of qubits. This is quite different from what happens in Shor's 
algorithm (compare with figure |5J|. 



necessary in this direction. In particular it should be very interesting to find out if exponentially large interference is 
a necessary condition for an exponential speed up of a quantum algorithm compared to its classical counterpart, or if 
the interference measure could be used to optimize existing algorithms or to conceive new ones. 
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